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Abstract 

We look at the covariant techniques and the ideas on constraints and gauge-invariance, which 
were recently employed in [1 to support earlier work by the same authors. That work was criticised 
in [2]. Using very simple and well known examples we show that, when adopted, the methods and 
views of [1] lead to basic-level mathematical problems, with analogous consequences for the physics. 
We provide a few simple rules that should help to avoid similar problems in the future. 

1 Introduction 

The subject of this note is related to the l+S-covariant study of the interaction between gravitational 
waves and large-scale magnetic fields in cosmology. Our aim is not to discuss this interaction again. 
The literature is available to anyone who might be interested. Here, we will look at the mathematics 
(as recently explained in [1]) on which the physics of [3] was based upon. 

The main issues have to do with the consistency of constraints, the use of 3-D commutation 
operators and gauge-invariance. The problems are technical but also elementary. The authors imposed 
inconsistent constraints, used incomplete commutation laws and presented as gauge-invariant results 
that are gauge-dependent by construction. We argue mainly by using simple and well known examples 
and provide a few basic rules that should help to avoid analogous problems in the future. Given that, 
the present brief communique may also have some pedagogical value. 



2 Constraints 



It was claimed (see § III.B in [3]) that, despite the presence of an inhomogeneous magnetic field, 
the electric component stays curl-free, provided it was so initially and the zero-order FRW model is 
spatially flat. The vanishing of curlE'^ does not depend on the scale, or the electrical conductivity, 
and applies to the second perturbative order. Even without looking at the mathematics, this claim 
sounds rather extreme. 

A constraint is consistent if, once imposed, it holds at all times with no need for further restrictions. 
To check whether curlE^^ vanishes, we must look for sources in the propagation equation (i.e. the first 
time-derivative) of curlEa- In [3j the authors considered instead the second time-derivative of this 
quantity, which on a spatially flat FRW background reads (see Eq. (29) in [3j) 



(curlEay = -I e(curlK)' + D^curl^ 
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curlE'a . 



Here, curlE'a = Sabc^'^E'^, Q is the background volume expansion, p and p are the matter density and 
pressure, A is the cosmological constant, is the 3-D covariant derivative and D-^ = D^D^ is the 
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associated Laplacian. The absence of explicit source terms in Eq. ([T]) misled the authors into claiming 
that curlE'a will remain zero, if it was so initially (see § III.B in [3j and also [Ij). This conclusion is 
incorrect, as the following simple example shows. 

Consider a perturbed FRW universe with dust. It is well known that the second time-derivative 
of the density perturbation (here represented by A) reads 

A = -^eA + ^pA. (2) 

Just like Eq. ([1]), the above contains no explicit source terms. Then, if we were to adopt [l^iSj, A will 
stay zero if it was initially zero and there will be no density perturbations0 We all know that this is 
not the case. To look for sources one must check the first time-derivative of A, namely the expression 

A = -Z. (3) 

This ensures that density perturbations are sourced by those in the expansion (represented here by Z). 
For the same reasons, when dealing with curli^^, one must look for sources in the first time-derivative 
of this quantity. The latter gives (see Eq. (1) in [2]) 

(curl^^))- = -ecurl^^) + TZ^^^ b\^-^ - D^S^) - curlj-i^) , (4) 

where TZab is the 3-Ricci tensor, Ja is the 3-current and the indices (1), (2) indicate the perturbative 
order of the variable. Thus, even if cvslEa = initially, it will not remain so once the magnetic and 
the current perturbations kick in (unless further constraints - on the sources - are imposed). The zero 
electric-curl claim made in [3] is unsustainable. 

Note that, in line with [3], the -B-field vanishes to zero order has a homogeneous first-order com- 
ponent (Bi^^^) and an inhomogeneous second-order part (^i^^). Thus, Ba = B^^^^ + B^ . Splitting 
the variables like that is largely avoided in covariant calculations, but we will adopt it here to agree 
with [1]. 



3 Commutators 

The 3-D commutators have the form 2DjqD{,] = DaD^ — D^Dq and monitor the commutation of 3-D 
covariant derivatives. How to use these operators became an issue, when pj sought to remove the 
source terms from the right-hand side of Eq. One can imagine situations where some of these 
sources are negligible relative to the rest, but not all three of them simultaneously. When the Laplacian 
is the only quantity left in the right-hand side of ([!]), in particular, we cannot ignore it (even on large 
- finite- scales) because it vanishes only asymptotically (i.e. at infinity). Here, we will focus on the 
curvature term because the argument used to remove this quantity suffers at the elementary level. 

The curvature term in Eq. ^ appears when analysing the quantity curlcurlSa by means of the 
covariant 3-D commutator Dj^Df,]. This operator applies to the full variable and the commutation 
takes place prior to any decomposition (into background, first-order perturbation, etc - the reasons 
will become clear below). It is argued in |l] (see Eq. (3) there), that we must instead split the variables 
first and commute their 3-D gradients afterwards. Splitting the magnetic field into Ba = B^^ + -Bi^\ 

(2) 

the authors involve only the inhomogeneous Ba -part and arrive at the second-order formula 

curlcurlBa = curlcurm^^) ^ _d2s(2) + TZfJ Bf^) = -B'^B^ , (5) 
^Initially, refers to the initial perturbed hypersurface and not to the FRW background, as [l] seems to suggest. 
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since TZ^^^ = (see Eq. (3) in [Ij). The term 'Tl^ah^\i) ® does not appear above, although it is also 
second order and it should have been included. 

We will accept, for argument's sake, the recipe leading to expression ([5]) and test it on a well known 
case. Assume a perturbed FRW universe containing a single perfect fluid with nonzero pressure. 
Following [1], we decompose the density a.s p = p^^^ + p^^\ where is the homogeneous (zero-order) 
part and p^^) the inhomogeneous perturbation. Then, in line with [Ij and Eq. ([5]), only p^^^ takes part 
in the operation and the linear commutator for the 3-gradients of p reads 

D[„D,]p = D[,D,]p(i) = -p(i)a;2) = , (6) 

since uj^^^ = 0. However, almost everyone using the covariant equations knows that the appropriate 
linear expression is 

D[,D,]p = D[,D,]p(i) = -p(°)a;i',) / , (7) 

obtained easily when we commute first and split afterwardsH The commutation recipe employed in |l| 
is incomplete. Note that if we were to use relation ([6]) instead of ([7]), the linear propagation equation 
of the vorticity would be given by 
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l^ab = -^Q^ab, (8) 

instead of the familiar expression 

t^ab = (^1 - ^ Swafe • (9) 

The extra term in the parentheses results from the right-hand side of Eq. ([7]) - is the square of the 
adiabatic sound-speed. Similar problems will emerge in several other expressions, if we were to follow 
the commutation rule of [Ij, and relativistic cosmology would require drastic revision. 



4 Gauge- invar iance 

An additional issue is that of gauge-invariance. In [H |3] the authors consider the aforementioned 
5-field, which has nonzero first-order value and is therefore gauge-dependent at second order by well 
known theorems. It is noted that the first-order field has 

5(1) + \ 65(1) ^ ^ (10) 

which makes the quantity 

Pa = Ba + \ QBa = B^^^ + ^ 65(2) ^ (H) 

a gauge-independent variable at second order. This is true, although there should be a 2G(i)5(i)/3 
term somewhere in the above (to simplify things lets say that 6(i) = 0). The real problem starts when 
the authors claim that they can integrate Eq. (jlip , with respect to Ba , and obtain gauge- independent 
results (see § IV in [Ij). This is an interesting idea, but it does not sound right. The gauge ambiguity 

(2) 

resides in Ba itself, by construction, due to the nonzero first-order value of 5^. 

Let us assume, for argument's sake again, that p^] are right and apply their principle to a familiar 
case. Consider a perturbed FRW universe containing dust with density p = pC^) + p(i). This quantity 
is the archetype of a gauge-dependent perturbation. To zero order, 

p(o)+0p(o) ^0, (12) 

^The reader might wonder why it should make a difference which operation (sphtting or commuting) is done first. It 
should not in principle, but if we first split and then commute the chances of doing something wrong increase considerably. 
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which means that the quantity 

£, = p + ep = p« + ep« (13) 

is a gauge-invariant first-order perturbation. Following [Ij, one can integrate the above and obtain 
gauge-independent results for p^^\ The same principle could be applied to every gauge-dependent 
perturbation. If this were true, there should be no need to introduce complicated variables to describe 
gauge-invariant perturbations and the gauge problem would have been a rather trivial one. Unfor- 
tunately, it is not that simple. Nevertheless, if the authors still believe in their case, they should 
communicate their findings without delay. 

5 Numerical results 

It is also worth mentioning the numerical evaluation of Eq. ()14p (see expression (7) in [T]). The latter 
gives the strength of a gravitationally amplified B-mode that passes through reheating and enters the 
horizon sometime in the radiation era 



B = Bn 



HJo \XhJrh\HJrh \^hJo\HJq\H 



RH 



(14) 



Here a/ H describes the shear anisotropy, and Xh are the scale of the i3-field and the horizon size, 
while the indices and RH indicate evaluation at the end of inflation and reheating respectively (with 
B = B^^^). The physics leading to the above result was explained in [2]. It was clarified there that the 
amplification is of the superadiabatic nature, takes place outside the horizon and defers from epoch to 
epoch. This view has also been adopted in [l], but the authors dropped the last term in the brackets 
claiming that it is of the third order. This is not so, as a bit of algebra can easily showH What is 
more important, is that the quantities in the brackets are numbers and what decides which stays or 
goes is their value. 

As explained in [2] (see Eq. (10) there), when there is substantial magnetic amplification during 
both the reheating and the radiation epochs, the third term in (I14p prevails. All these are of uncertain 
practical value, however, as long as the gauge-issue remains unresolved. 

6 Conclusions 

We close by summarising this communique in four basic rules: (i) To check the consistency of a 
constraint we look at the first (not the second) time-derivative of the constrained quantity; (ii) To 
avoid mistakes, we first commute the 3-D gradients and then split the involved variables (not the other 
way around); (iii) The gauge-ambiguity and the arbitrariness are in the variable that describes the 
perturbation; (iv) When evaluating sums like that of Eq. (fT4|) . one should always ensure that all the 
terms are properly accounted for. 
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^It helps to recall that the amplified B-field at the end of reheating, acts as the background source for the gravito- 
magnetic interaction in the radiation era (see [2]). 
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